function OpenCRTBP_u(u)
%% Constants
figure
hold on
[L1,L2,L3,L4,L5] = librationPoints(u);
Lse = [L1,L2,L3,L4,L5];
for ii = 1:5
    plot3(Lse(1,ii),Lse(2,ii),Lse(3,ii),'bo','MarkerFaceColor','b','MarkerSize',2)
end

plot3(1-u,0,0,'rx');
plot3(-u,0,0,'mx');
xlim([-1.5,1.5]);
ylim([-1.5,1.5]);
axis square
xlabel('X [DU]')
ylabel('Y [DU]')
zlabel('Z [DU]')
set(gca,'FontSize',20) 
end